Pulsar motions from neutrino oscillations induced by a violation of the equivalence 

principle 



M. Barkovicht, H. Casini*, J.C. D'OhW, R.Montemayor* 
1 Institute: de Ciencias Nucleates, UNAM, Ap. Postal 70-543, 04510 

Mexico DF, Mexico 
* Theoretical Physics, University of Oxford, 1 Keble Road, Oxford 
OX13NP, United Kingdom 
* Instituto Balseiro and CAB, Universidad Nacional de Cuyo and CNEA, 
8400 Bariloche, Argentina 

We analize a possible explanation of the pulsar motions in terms of resonant neutrino transitions 
induced by a violation of the equivalence principle (VEP). Our approach, based on a parametrized 
post-Newtonian (PPN) expansion, shows that VEP effects give rise to highly directional contribu- 
t-H ' tions to the neutrino oscillation length. These terms induce anisotropies in the linear and angular 

momentum of the emitted neutrinos, which can account for both the observed translational and rota- 
tional pulsar motions. The violation needed to produce the actual motions is completely compatible 
with the existing bounds. 
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It is very difficult to obtain precise evidence on the characteristics of the gravitational interaction beyond the range 
where the Newtonian approximation holds. Only systems with very large densities of mass in rapid motion can provide 
suitable laboratories for such a phenomenology. One well known example is the orbital behavior of binary pulsars, 
■ which gives support to the production of gravitational waves. Type II supernovas are another interesting scenario. In 
f^- \ this case the intense neutrino flux produced during the gravitational collapse can be sensible to subtle characteristics 
. of the gravitational interaction. In this letter we analyze some effects on this flux that could test a possible violation 
' of the equivalence principle. 

Perhaps one of the most intringuing characteristic of the pulsar dynamics related with the supernova stage is their 
anomalous proper motions. There is strong observational evidence that translational velocities of pulsars include a 
significant component from kicks given when they are formed jlj] . Several mechanisms have been proposed to explain 
such kicks, but none of them is completely satisfactory ||,||||. Recently it has also been pointed out that the observed 
rotation periods are several orders of magnitude shorter than the predictions for the cores of the protoneutron stars 
§. Thus the spin of the pulsars are probably produced by the same mechanism that gives them their translational 
velocities during the formation stage. Moreover, there is significant observational evidence that seems to indicate a 
polarization of the motion of young pulsars along a direction near the plane of the galaxy Q . This correlation could 
mean that kicks involve a characteristic length at least of the order of the galaxy radius, which is very difficult to 
. . , explain on the basis of the proposed mechanisms. 

An appealing possibility that could account for the translational kick is a 1% anisotropy in the momentum carried 
by the neutrinos emitted during pulsar formation. However small, such anisotropy is not easy to obtain. Kuscnko and 
Segre (KS) || have proposed a mechanism based on the deformation of the resonance surface when neutrinos undergo 
matter oscillations in the presence of a magnetic field J8j . Unfortunately the necessary magnetic field is relatively high, 
B > 10 15 G Furthermore, the condition that the resonance surface has to lie between the neutrinoespheres 

implies m v ~ 100 eV. The existence of such heavy neutrinos is cosmologically ruled out unless they are unstable. 

A less orthodox mechanism for neutrino oscillation was proposed several years ago. It requires a flavor dependent 
coupling of neutrinos to gravity |pd| , and no neutrino mass. Consequences of such a violation of the equivalence 
principle (VEP) in the neutrino sector have been analyzed in a number of papers @. In particular, in Ref. [§ it 
was applied to the problem of the translational motion of pulsars . In this case the desired kick can be achieved with 
massless neutrinos, but the intensity of the magnetic field is similar to the one required by the KS mechanism. 

In this work we propose a purely gravitational explanation for both the translational and rotational motion of 
pulsars, where the neutrino oscillation and the momentum anisotropy are induced by VEP effects and that do not 
rely on the magnetic field of the protostar. We work within the framework of a generalized parametrized post- 
Newtonian (PPN) formalism previously applied to the solar neutrino problem |T4| ], that naturally includes the 
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effect of a preferred reference system. Our approach generalizes the usual VEP scheme |12 , by including the effect 
of potentials of the next PPN order to the Newtonian potential U, and a tensorial potential of the same order than 
U. In principle all these terms should be present if the equivalence principle is violated. In this context the neutrino 
oscillations are a manifestation of a VEP effect, and the momentum anisotropy is signature of the preferred reference 
system. The accuracy of the equivalence principle may be characterized by limits in the differences of the PPN 
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parameters for different neutrinos. As we show, violations of the equivalence principle consistent with the present 
bounds generate the necessary kicks to produce the observed pulsar motions. 

The linearized Dirac equation for massless neutrinos in a static gravitational field leads to the dispersion relation 
0: 



1 + Pi ~ TjhtjPipj - ^h 00 U 



(1) 



where the W v fields are defined by = rf ,v + h^ v , referred to the Minkowskian metric. In deriving this relation we 
have neglected the spatial derivatives of the gravitational potentials, which is justified for neutrinos in astrophysical 
systems. Up to third order in the velocity of the source w we can write (G — h = c = 1): 



h 00 = 2j'U + 0(w 4 ), 

7 1 1 

Ki = - ^ AiV< - -A 2 Wi + (a 2 - -a^v.U - a 2 v J U jl 

hij = 2 7 C% + TUn +0(w 4 ). 
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The adimensional parameters of the PPN expansion are 7, 7', Ai, A2, T, v , ai, and ol%. The parameters ol\ and 
a 2 vanish in Lorentz covariant theories, but if there exists a preferred reference frame, characterized by a velocity v, 
they should be non null. 

The general expressions for the potentials U, Vi, Wi and Uji can be found in Ref. |Q. In the present case, the 
source of the gravitational field is the protoneutron star. Considering a spherical configuration and a rigid rotation, 
the PPN potentials become 



U = 4tt 



dr r 
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Or - r') + —6 (r' - r) 

V*' 



P{r') , 



Uij = f l f :j I(r) + 5ij J(r) , 
V, = Wi = Wl J{r) . 

where p(r) is the mass distribution of the star and Wi — eijk^ljTk- Here O is the angular velocity and 
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In presence of VEP all the PPN parameters can depend on the flavor numbers. We assume that deviations from 
a metric theory are small, so that in a very good approximation there is a common coordinate frame for all flavors. 
Since the parameters are flavor dependent, distinct neutrinos will undergo different phase shifts when passing through 
the same sector of space. In the presence of neutrino mixing phase shift differences become observable as neutrino 
oscillations. For simplicity, in what follows we consider two neutrino flavors, v e and or v r . They are supposed to 
be linear superpositions of the gravitational eigenstates v\ and v 2 , with a mixing angle 9 g . Along the neutrino path 
flavor evolution is governed by 



Ao 
2 



cos26> 9 sin 29 g 
sin 29 g cos29 g 



dr V^m 

with Ao = E 2 — E . For a rotating protoneutron star we have 

A = |-(<5 7 ' + 5j)U - 5TJ - 6TI(r ■ p) 



(10) 
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(Sa 2 — -8a,\)U — 5a 2 J 



--(7<5Ai +SA 2 ) ,m xr-p}E, 



v • p — 5a 2 I(r ■ v)f ■ p 



(11) 



where E = p is the neutrino energy, 6j = j 2 — 7 1 , and the same for the difference between the other PPN parameters. 
Here Aq plays the same role as the quantity (m| — m\)/2E in the mass mechanism for neutrino oscillations. Note 
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that in our case the potentials depend on r and hence Ao = Ao(r). Terms with v appear whenever a preferred 
frame exists. In principle v could also depend on the gravitational flavor, but the observed position offset for pulsar- 
supernova remnant pairs can be interpreted as the existence of a translational effect associated to a preferred 
direction. For this reason we take v as a flavor independent parameter. Its action is analogous to the one produced 
by a magnetic field in the KS mechanism. 

As is well known, neutrino oscillations in matter differ from the oscillations in vacuum. The interaction of neutrinos 
with the background modifies their dispersion relations, and under favorable conditions leads to the MSW phenomena 
of resonant flavor transformation. If electrons are the only leptons present in the medium, the term ^=iVe(r)<73 has 

to be added to the matrix in Eq.(|ni|), where 03 is the Pauli matrix and N e (r) denotes the electron number density. 
The resulting Hamiltonian can be diagonalized at every point by a local rotation, with the mixing angle in matter 
9 m (r) given by sin26> m (r) = A ^ sin26> g , with 



A(r) = J (A (r) cos20 g - V2G F N e (r)) + (A (r) sin20 ff ) 



(12) 



There is a resonance when the diagonal elements of the Hamiltonian vanish, i.e. when \/2GFN e (rn) = 
Ao(r_R,) cos26 g . The efficiency of the flavor transformation depends on the adiabaticity of the process, which is 
characterized by the parameter 
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where the scale heights are h 
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lniV e and h~^~ = ^ In Ao. The transition will be adiabatic whenever k 1. 
The translational kick comes from the anisotropy in the radial momentum carried by neutrinos emerging from the 
resonance surface. The resulting effect on the motion of the pulsar is obtained integrating over all the surface. To 
this integration will only contribute the radial component of p. Therefore, to estimate the translational kick we use 
a simplified situation with a purely radial neutrino flux. In this case, 



A = [A(r) + B(r)vcosx\E , 
where x is the angle between r and v. The functions A(r) and B(r) are given by 

A = -(<V + 5~/)U -5T(I + J) , 
B = (Sa 2 - -Sai)U - 5a 2 (I + J) . 

The radius of a point on the distorted resonance surface can be written as tr = r Q 
of the unperturbed resonance sphere r is determined by 



A(r ) = 



V2G F N e (r ) 



cos 28 n 



E 
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S cosx (S <C r a ). The radius 
(17) 



and 
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(18) 



where we keep only the terms linear in 8, and — — In A(r). 

At the moment there is no agreement about the details of the production of a kick by a distorted neutrinosphere. 
To explore the possibilities of the VEP mechanism we will now consider this effect in the context of the main 
neutrinosphere models proposed. 

For a hard neutrinosphere model in thermal equilibrium as considered in Refs. (3|4|U, the momentum asymmetry 
in the v direction is generated by the emission at points with different temperatures on the resonance surface 
Ap/p 

,1/3 -' :x -|T/x 



^hj}5, where h T l = g~lnT. In the case of a quasi-degenerate gas of relativistic electrons with a constant 



chemical potential fi e w (37r 2 7V e 
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with Q = ^fp?, where 7/ = n e /T is the degeneracy parameter for the electrons and A = hji/(hA — hjy e ). Another 
possibility is to assume that the electron fraction Y e remains constant and p ~ T 3 |9). In this case hjsr e ~ and 
Q = ^A. 

A different kick model in the literature uses a soft neutrinosphere [§,[l0| . In such a case there is an important 
reduction in the anisotropy given by the ratio p Q j p c of the density at the resonance and the density at the core. 
The momentum asymmetry can also be written as in Eq. (|l9|) , with Q — p h^ c A/18m c , where m c — f rs pdr is the 
integral of the mass density between the central core and the surface of the star. In all the cases considered above the 
adimcnsional parameter Q depends only on the specific model and the remaining factors contain the PPN parameters. 

For a quantitative estimation of the effects of VEP on the neutrinosphere, we use the density profile p(r) = p c for 
r < r c and p(r) — p c (r c /r) n for r > r c ||. We take p c = 8 x 10 14 g/cm 3 , r c = 10 km, and 5 < n < 7 that give a 
good description of the supernova SN1987A |jq| . The resonance surface has to lie below the v e neutrinosphere and 

above the neutrinosphere. If we take p Q ~ 10~ n (//cm 3 and Y e ~ 0.1, then we obtain [S^ + S^ + 0.95 ST\ cos 28 



—6 x 10~ 10 . For ST = our result agrees with the one obtained in Ref. ||. As pointed out in this work the adiabaticity 
condition is achieved provided that 8 g > 10~ 4 , h^ 1 < ^aJ' anc ^ hence A ~ 1 for every value of n. The value of the 
momentum asymmetry is 

Ap 

P 



-Q(5ai -O.l5a 2 )vcos26 g x 10 9 . (20) 



For T — 3MeV, Q ~ 0.1 in the hard neutrinosphere models, and Q ~ 4 x 10 -5 for a soft neutrinosphere. Taking 
Sai ~ Sa2 ~ (5a, and requiring Ap/p ~ 0.01, we obtain z;i5a ~ 10 -10 and w<5a ~ 10 -7 , respectively. 

We now analyze the effect of the non radial component of the neutrino momentum. When S~2 = 0, at a given point 
of the resonance surface the emitted neutrinos have an azimuthal symmetry respect to the position vector. For non 



vanishing angular velocity of the protoneutron star, the last term in Eq. (11) brakes this symmetry and produces an 
angular acceleration of the star. To make a perturbative estimation of this effect we ignore the dependence of A D on 
v and adopt a very simple model of a hard resonance surface at tq + Sr. From the resonant condition we get 



C 

Sr = — Aft,jv c fi'ixp 



(21) 



where C(r) = — | (7<5Ai + SA 2 ) J (r). 

Neutrinos emitted in different directions come from regions at different r and therefore at different energies. Hence 
they have different angular momenta. If we adopt the Stefan-Boltzmann law for the neutrino flux at the resonance 
surface, a neutrino emitted in a direction p has a momentum p = E D (1 + Ah^Sr), where E a = E(r a ). Therefore it 
carries an angular momentum 

1 =r E {r x p) [l + Ah^Sr] . (22) 

By integrating at each point of the resonance surface over all directions and also over all the points, we compute the 
angular momentum gained by the star. Because of the symmetry of the system the resulting angular acceleration 
points along the rotational axis. The time derivative of the total angular momentum can be expressed as 

• CA h N ■ lo de sin 9 io f d6 ' $T V sin e ' r o (« x f ■ p) 

L = -J±-£nr ; i '—, (23) 

3ttA h T ^ dO sin 8 / * d6> sin 8' 

where £ is the energy carried by the neutrinos per time unit, and a factor g has been included to take into account 
that although six neutrino and antincutrino species are radiated, only one comes from the distorted neutrinosphere. 
In the latter expression 8 is the angle between the radius vector f and the angular velocity fl, while 8' and ip' are the 
spherical coordinates for p taking f as the z axis. From Eq. (EF 



ne>-«.«p(£jf^s*), (24) 

where / and f2 are the momentum of inertia and the initial angular velocity of the protostar. It should be noted that 
the rotational kick does not require a velocity v associated to a preferred frame. 
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As an example, let us consider the density profile introduced above. Assuming that all the quantities in the integrand 
except £ are constant during the cooling period and taking A£ ~ 3 x 10 53 erg, the angular velocity after the angular 
kick is 



Jl/ ~ n o exp 



^UA 1 + W)^10* 



(25) 



where O.f < £ < f for 5 < n < 6. The ratio -^f- is a model dependent quantity of the order of unity. Values usually 

considered are h^ e < hx/3. The star angular velocity will increase or decrease depending on the sign of SAi + jSA2- 
If we accept that typical initial angular velocities of the protostar are Sl ~ 0.010/, the VEP parameters must be in 
the range 10 -6 < SAi + jSA2 ^ 10 -8 to reproduce the observed values for SI/. 

To estimate the order of magnitude of the translational and rotational accelerations, we have assumed the corre- 
sponding kicks decoupled one from the other. In a more realistic situation the rotational motion could produce an 
average in the translational kick. This effect depends on the relation between the characteristic time of reaccommo- 
dation of the neutrinosphere and the period of rotation of the star. The anisotropy axis coincides at every time with 
v and is not affected by the rotation, but the temperature of the resonance surface could change. In a soft neutri- 
nosphere the deformed resonance surface changes the atmosphere opacity over the core and induces a temperature 
anisotropy in the core-atmosphere interface, which in turn affects the neutrino flux. As the star rotates the resonance 
surface also rotates with respect to the rest frame of the star inducing a time-changing opacity over the core region. 
Therefore we have here to consider the characteristic thermal response time of the system, which is of the order of 
a few hundred miliseconds ]To[ | , in contrast with the pulsar period. Thus, in this case we can expect an averaging of 
the translational kick. This effect tends to cancel the component orthogonal to the rotational axis and develops a 
correlation between the translational kick and the axis of rotation. In the case of a hard neutrinosphere the energy 
flux depends on the temperature at the point from which the neutrinos are radiated from the surface of resonance. 
The atmosphere here has enough heat capacity to act as a thermal reservoir with a radius dependent temperature. 
Therefore, there is no effective average and there is no correlation between the translational kick and the rotational 
axis. 

For simplicity, we have assumed that the only mechanism responsible for the pulsar motion is VEP. If this were the 
only cause for the translational velocity, then all pulsar velocities should show a certain correlation driven by the v 
parameter. This correlation will be more or less accentuated depending on how hard or soft the neutrinospheres are, 
and also could be blurred by the presence of other kick mechanisms besides the one here considered. 

In conclusion, we have shown that resonant VEP neutrino oscillations may be responsible for both the translational 

and rotational motion of pulsars. Since this mechanism works even for massless neutrinos, it does not clash with 

cosmological bounds. The strictest boundaries known at present in the neutrino sector are given by accelerator 

experiments, mainly from CCFR, which correspond to the highest tested energies fllTfl . These experiments are sensitive 

to large mixing angles, because they have no access to the MSW effect. The exclusion region for these experiments 

extends down to sin 2 29 > 2.10 -3 for Venu^ and sin 2 29 > 0.2 for v e nu T , independently of the value of Ao. Therefore, 

the parameter region relevant for the neutrino resonance in neutron stars, taking 10~ 4 < 9 g < 10~ 3 for v e nu^ or 
i-4 ^ a ^ m-i 



10 < 9 g < 10 for v e nu T , is well outside the range tested by accelerators Jl8|]. With respect to atmospheric 
neutrinos, they are not affected by these small mixing angle oscillations, and the MSW effect in the solar neutrinos 
corresponds to a medium of much lesser density, and thus the involved parameter sector is very different [p"i|| . In this 
way the kick pulsar physics gives access to a new phenomenological sector of VEP effects. 
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